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I.  INTRODUCTION 

The  purpose  of  Phase-I  of  the  En  Route  Radar  Weather  Extraction  Program 
was  to  develop  and  recommend  techniques  needed  to  generate  accurate  weather 
level  signals  In  the  en  route  environment.  Bap  has  Is  was  to  be  placed  on  the 
accurate  measurement  of  precipitation  reflectivity  in  the  presence  of  ground 
clutter.  The  results  of  previous  work  in  this  area  by  the  FAA,  National 
Severe  Storms  Laboratory  (NSSL) ,  FAA  Technical  Center  (FAATC) ,  and  the  Johns 
Hopkins  Applied  Physics  Laboratory  were  to  be  used  as  technical  background  for 
the  study.  The  techniques  used  by  the  FAA  in  Ref.  1  and  by  the  National 
Weather  Service  in  Ref.  2  to  provide  six  levels  of  weather  information  were 
investigated  as  a  baseline  for  the  feasibility  and  utility  of  this  approach. 

It  should  be  emphasized  that  the  principal  objective  of  the  Riase  I 
studies  was  to  examine  a  variety  of  candidate  techniques  using,  when  possible, 
actual  Moving  Target  Detector  (MTD)  data  to  assess  the  performance  limitations 
due  to  ground  clutter.  The  resulting  candidate  weather  radar  front  end  and 
signal  processor  has  been  based  on  this  practical  experience,  but  has  not  as 
yet  been  fully  implemented  and  evaluated  in  an  operational  environment. 

Currently,  weather  information  is  obtained  from  the  en  route  radars  and 
sent  to  the  Air  Route  Traffic  Control  Center  (ARTCC)  via  the  weather  and  fixed 
map  unit  (WFMJ)  and  common  digitizer  (CD).  The  WFMU  thresholds  the  radar 
video  and  the  CD  generates  digital  messages  giving  the  azimuth  and  range  of 
weather  threshold  crossings  and  transmits  this  information  to  the  ARTCC  where 
it  is  used  to  produce  a  digital  display  of  the  weather.  Die  measurement  of 
this  weather  data  suffers  from  inaccuracies  due  to  circular  polarization  (CP) 
attenuation,  sensitivity  time  control  (STC)  attenuation,  and  the  velocity 
response  of  the  moving  target  indicator  (MTI)  circuitry  used  to  eliminate 
ground  clutter. 

Modern  radar  digital  signal  processors  have  been  developed  which  tend  to 
eliminate  the  effects  of  ground  clutter,  angel  false  alarms,  and  precipitation 
echoes  in  order  to  furnish  displayable  target  reports  of  moving  aircraft 
only.  These  processors  are  suitable  for  use  with  the  FAA  L-Band  radars  used 
to  perform  surveillance  functions  in  en  route  air  space.  To  eliminate  false 
alarms  due  to  precipitation  echoes,  adaptive,  linear  constant  false  alarm  rate 
(CFAR)  techniques  have  been  used  in  the  Lincoln  MTD  radar  signal  processors. 
Under  this  program  it  was  planned  to  use  real-time  MTD  measurements  in  range- 
angle-velocity-space  after  normalization,  to  estimate  the  reflectivity  factor 
(i.e.,  z-value)  of  precipitation  returns. 

The  technique  developed  for  use  with  the  en  route  radar  equipped  with  MTD 
processing  uses  a  two-level  ground  clutter  map  to  aid  the  weather  extraction 
processf^J.  Basically  the  signal  processor  generates  two  sets  of  weather 
threshold  crossings.  One  uses  the  total  received  energy  of  the  radar  returns 
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and  the  other  uses  only  a  linear  combination  of  the  non-zero  doppler  filter 
outputs.  For  high  clutter  cells,  the  non-zero  frequency  filter  outputs  are 
used.  Otherwise  the  total  received  energy  is  used.  These  data  are  then  inte¬ 
grated  over  three  scans  and  spatially  smoothed  to  generate  a  contoured  weather 
map.  This  process  is  utilized  for  two  selectable  z-value  levels  of  precipita¬ 
tion  reflectivity. 

II.  PROGRAM  TO  DATE 

A  program  to  enhance  the  clutter  performance  of  terminal  (ASR)  and 
en  route  (ARSR)  radar  systems  by  introducing  MTD  processing  has  been  underway 
for  some  time  (Ref.  3).  As  a  by-product  these  MTD  processors  were  adapted  to 
produce  two  calibrated,  contoured  levels  of  rain  reflectivity  within  coverage 
of  the  sensor.  It  was  believed  that  a  six-level  weather  extraction  processor 
might  reside  in  the  back-up  channel  MTD  processor,  and  use  the  fine-grain, 
temporally  smoothed  ground  clutter  map  (maintained  by  the  MTD  for  zero-veloc¬ 
ity  target  thresholding  against  clutter)  to  subtract  from  the  current  measure¬ 
ment  of  weather  plus  ground  clutter,  producing  a  cell-by-cell  estimate  of 
weather  only.  Early  efforts  focussed  on  exploring  this  possibility  by  devel¬ 
oping  smoothed  ground-clutter  maps  using  the  ASR-7/MTD  at  Lexington,  MA;  the 
ASR- 7/ MTD  at  Burlington,  VT  (BTV) ;  the  FPS-67/MTD  at  Bedford,  VA  (BVA);  and 
the  instrumented  ASR-8  at  the  FAATC,  Atlantic  City,  N.J.  This  study  was  con¬ 
ducted  during  the  period  January-September,  1979.  However,  by  the  fall  of 
that  year  we  were  convinced  that  this  technique  could  not  be  made  to  support 
weather  extraction  satisfactorily  for  the  following  reasons: 

1)  The  estimate  of  ground  clutter  exhibited  higher  than  expected  variance 
because  the  "large"  cross  section  cells  consisted  of  only  a  few  singularities, 
and  the  single-frequency  radar  measurement  coupled  with  this  feature  produced 
low  frequency  variations  which  were  tracked  by  the  single-pole  low-pass  filter 
used  for  smoothing.  The  variations  in  mean  radar  cross  section  were  as  great 
as  +6  dB,  and  it  was  determined  that  rejecting  the  ground  clutter  return  at 
zero-velocity  would  be  required  in  order  to  develop  an  estimate  of  the 
non-zero  velocity  components  of  the  weather  in  the  cell(s)  of  interest.  An 
example  of  ground  clutter  amplitude  variation  with  the  BVA  FPS-67  antenna 
stopped  is  shown  in  Fig.  II-l. 

2)  Those  cells  that  did  not  exhibit  the  above  variations  (wooded  areas) 
did  exhibit  very  slow  variations  of  up  to  6  dB  from  day  to  day,  seemingly  as  a 
function  of  wet/dry  conditions.  The  MTD  processed  approximately  5x10^ 
range/azimuth  clutter  map  cells  (1/16  nmi  x  0.7°),  and  although  only  10Z 
(above  +30  dB  were  affected  by  the  above  problems,  this  was  con¬ 
sidered  to  be  sufficiently  severe  (statistically)  to  prevent  recommending  that 
this  technique  be  Implemented. 
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zero-velocity  filter  output  with  antenna  stationary 


3)  Use  of  the  ARSR/MTD  processor  for  surveillance  and  multi-level  weather 
extraction  processing  complicates  the  antenna  polarization,  range,  and  radar- 
STC  normalization  of  the  weather  thresholds.  If  the  system  were  to  be  hard¬ 
wired,  and  pre-computed  normalized  weather  thresholds  implemented  in  ordered 
range,  then  radar/weather  processing  could  co-exist.  However,  the  radar  MID 
sub-system  was  to  be  implemented  using  fault-tolerant  architecture  which 
allowed  for  restructuring  of  the  processor  to  nonrrange  order  under  software 
control.  Thus  to  maintain  this  fault-tolerant  architecture  would  require  on- 
the-fly  re-computation  of  weather  thresholds  module-by -module. 

4)  Finally,  use  of  the  Doppler  filter  bank  optimized  for  aircraft  detec¬ 
tion  for  weather  extraction  processing  was  deemed  inappropriate.  This  is  true 
because  a  large  fraction  of  the  precipitation  cells  will  have  a  velocity  spec¬ 
trum  centered  at  zero  radial  velocity.  Hence  it  would  be  necessary  to  ex¬ 
clude  a  large  velocity  interval  around  zero-velocity  to  minimize  false  air¬ 
craft  target  declarations.  This  low  frequency  filtering  would  suppress 
weather  returns  (along  with  the  ground-clutter  returns)  and  result  in  distor¬ 
ted  weather  contours. 

The  above  considerations  led  us  to  investigate  the  possibility  of  imple¬ 
menting  a  separate  weather  channel  processing  system  which  could  operate  in 
parallel  with,  and  independently  of  the  MTD  system.  This,  of  course,  gives  up 
access  to  the  MTD  zero-velocity  clutter  map,  but  this  was  no  longer  considered 
to  be  useful  for  extracting  weather  attributes  from  weather-plus-ground-clut¬ 
ter  radar  cells. 

III.  CANDIDATE  RADAR  FRONT-END  AND  SIGNAL  PROCESSOR  FOR  WEATHER  DATA 
EXTRACTION 

The  remainder  of  this  working  paper  deals  with  a  candidate  weather  ex¬ 
traction  processor  for  use  with  either  the  ASR-MTD  or  ARSR-MTD  radar  systems 
being  developed  for  FAA-ATC  aircraft  surveillance.  Also,  Included  in  the 
Appendices  are  a  comparison  of  the  weather  return/receiver  noise  ratios  for 
three  candidate  radars,  data  to  place  clutter  and  precipitation  levels  in 
perspective,  comments  on  the  operation  of  the  ARSR/MTD  in  the  field,  and  data 
on  candidate  clutter  filters. 

A.  Front-End 


A  number  of  front-end  configurations  were  considered,  during  the  course 
of  this  study,  and  all  could  be  made  to  support  weather  extraction.  The 
candidate  front-end  configuration  shown  in  Fig.  III-l  is  the  preferred 
arrangement,  as  it  decouples  the  weather  channel  from  gain  variations  induced 
by  changes  in  radar-STC  curves.  It  does,  of  course,  introduce  a  1  dB  loss  in 
the  primary  channel,  but  this  is  considered  acceptable  as  tilting  the  antenna 
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down  will  more  than  compensate  for  this  loss  In  the  MTD  channel.  Translation 
to  IF,  quadrature  video  processing,  and  A/D  converter  sampling  of  I  and  Q  var¬ 
iables  is  the  same  as  specified  for  the  ARSK-MTD.  Die  weather  channel  STC 
should  be  operated  at  a  l/R^  rate,  to  normalize  all  weather  signals  as  shown 
in  Fig.  II1-2. 

B.  Digital  Signal  Processor 

Die  high-speed  digital  signal  processor  (ground  clutter  filters)  performs 
a  function  similar  to  the  MTD  Doppler  filter  bank,  but  is  tailored  to  the 
weather/ ground  clutter  problem,  and  operates  slightly  differently.  A  major 
element  of  this  study  was  to  characterize  the  ground  clutter  statistics  as 
they  applied  to  the  weather  extraction  problem,  and  to  develop  a  set  of 
reasonable  filters  which  would  permit  accurate  estimation  of  weather  in  the 
presence  of  ground  clutter.  One  obvious  choice  was  an  all-pass  filter  which 
would  be  optimum  for  large  Wx/Ci  range-azimuth  cells,  since  this  filter  intro¬ 
duces  zero-error  into  the  weather  estimate,  and  one  of  the  two  MTD  filters  for 
weather  was  all-pass.  Most  cells,  for  most  weather  levels,  will  be  processed 
using  this  filter.  Possibly  10Z  of  the  cells/levels  will  require  high-pass 
filters  to  improve  the  Hx/Ci  ratio  prior  to  estimating  the  precipitation 
reflectivity.  Another  simple  filter,  which  might  be  the  next  choice,  is  a 
four  pulse  Finite  Impulse  Response  (FIR)  filter.  This  filter  would  operate  on 
the  8  pulses  in  a  CPI  (I  and  Q  samples  for  each  pulse)  to  yield  A  values  of 
steady  state  output.  The  weather  level  would  be  estimated  from  these  A 
values.  The  transfer  function  of  this  filter  is  shown  in  Fig.  III-3.  The 
clutter  rejection  for  this  filter  is  approximately  -12  dB,  which  would  handle 
the  clutter  at  ranges  greater  than  20  tan  for  most  weather  levels  of  interest. 
This  filter  will  reject  some  of  weather  spectra  which  fall  in  the  range  from 
-2  m/s  to  +2  m/s  and  hence  should  only  be  used  when  necessary.  The  small  per¬ 
centage  of  range-azimuth  cells  which  require  more  than  12  dB  of  clutter  rejec¬ 
tion  could  use  2-pulse  or  3-pulse  cancellers  such  as  are  described  in  Appendix 
B. 


During  the  investigation  of  candidate  filters,  the  clutter  map,  as  re¬ 
corded  at  BVA,  was  run  against  four  filters,  and  a  conservative  clutter  rejec¬ 
tion  requirement  used  to  examine  the  resulting  weather  maps.  The  plots  of 
Figs.  III-A  to  -11  show  the  location  of  cells  requiring  the  use  of  a  given 
level  of  filtering  to  achieve  detection  of  weather  signals  of  a  specified  in¬ 
tensity.  Die  maps  were  prepared  using  alternate  CPIs,  yielding  256  azimuths, 
and  the  limiting  range  of  clutter  at  BVA  of  approximately  80  nml.  There  were 
A0960  range-azimuth  cells  (l.A°  x  1/2  nmi  x  80  nmi)  processed  in  covering  the 
weather  (636  one-eighth  nml  gates  were  processed  in  one  nml  blocks),  and  a 
choice  of  four  clutter  filters  was  available  for  each  cell.  The  clutter  can¬ 
cellation  of  the  filter(8)  used  in  each  mapping  exercise  is  given  in  Table 
III— 1 ;  the  response  of  each  filter  is  as  described  in  Appendix  B. 
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Fig.  III-3.  Frequency  response  of  candidate  FIR  clutter  filter 


20  NMI  RANGE  RING8 

Fig.  III-4.  57  dBz  cells  requiring  filter  #1. 
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Fig.  III-7.  46  dBz  cells  requiring  filter  #3. 
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Fig.  III-8.  46  dBz  cells  requiring  filter  #1. 
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Fig.  III-9.  41  dBz  cells  requiring  either  filter  #2 
or  filter  #3,  mostly  #2. 
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Fig.  III-ll.  30  dBz  cells  requiring  filter  #3. 


TABLE  III-l 

FATHER  MAPPING  EXERCISE  SUMMARY 


Precipitation  Level 

Cells  Requiring  Filter 

Figure 

dBz 

Descriptor 

III- 

57 

Intense 

1  (12  dB) 

4 

50 

Very  Strong 

2  (20  dB) 

5 

50 

Very  Strong 

1  (12  dB) 

6 

46 

Strong 

3  (40  dB) 

7 

46 

Strong 

l  (12  dB) 

8 

41 

Moderate 

2  or  3  (20  or  40  dB) 

9 

41 

Moderate 

3  (40  dB) 

10 

30 

Weak 

3  (40  dB) 

11 

The  analysis  was  carried  out  using  rain  intensity  values  of  57,  50,  41, 
and  30  dBz  corresponding  to  a  range  from  intense  to  weak  precipitation  (Table 
III-l  summarizes  the  mapping  exercise  where  the  NWS  levels  of  precipitation 
are  as  defined  in  Fig.  A-6) .  For  the  Intense  precipitation,  most  of  the  cells 
can  be  processed  by  the  all-pass  or  Type  0  filter;  Fig.  III-4  shows  the  few 
cells  requiring  Filter  Type  1,  a  four  pulse  canceller.  With  less  conservative 
criteria  even  fewer  cells  would  require  this  filter.  The  57  dBz  intensity 
signals  did  not  require  the  use  of  the  20  or  30  dB  filters. 

As  the  rain  signal  decreases  in  intensity,  more  clutter  filtering  is 
required  as  illustrated  in  Fig.  I1I-5  and  -6  which  show  the  cells  requiring 
the  use  of  the  20  and  12  dB  filters  respectively.  There  were  few  cells  which 
required  the  20  dB  filter.  Further  decreasing  the  Intensity  to  46  dBz 
requires  more  robust  filtering  as  shown  in  Fig.  III-7  and  -8,  where  both  the 
46  dB  filter  is  engaged  for  a  few  cells,  and  a  large  number  of  cells  require 
the  12  dB  filter.  The  20  dB  filter  is  required  for  processing  the  41  dBz 
intensity  level  as  is  the  40  dB  filter.  See  Fig.  III-9  and  -10.  Again  the 
number  of  cells  requiring  the  large  clutter  attenuation  is  small.  Finally,  at 
the  weakest  intensity  level,  30  dBz,  the  40  dB  filter  is  required  to  reject 
clutter  in  a  large  number  of  cells. 

It  is  believed  that  a  slight  change  in  the  clutter  rejection  criteria 
could  reduce  the  need  for  the  high-attenuation-value  filters.  Parametric 
analyses  were  not  carried  out  for  a  number  of  other  clutter  rejection  options. 

Each  filter  operates  on  all  data  through  weather  level  thresholding,  at 
which  point  the  site-dependent  map  selects  one  of  four  threshold  declarations 
for  each  1/2  nmi  range  Interval.  In  the  absence  of  rain,  weather  thresholds 
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are  not  declared,  even  though  random  threshold  crossings  occur.  Ihe  map  is 
created  in  clear  weather,  and  iterated  to  reduce  false  declarations  to  "near 
zero". 

Scan-to-scan  smoothing,  and  contour  s tart/stop  detection  is  as  specified 
for  the  ASR-9/MTD  two-level  weather  extraction  processor.  Nominally,  all  five 
contour  maps  could  be  continually  maintained  and  updated,  but  it  is  likely 
that  only  two  would  be  displayed  simultaneously,  and  it  would,  of  course,  be 
more  cost  effective  to  develop  only  the  contours  which  are  to  be  displayed. 
Still,  this  represents  straightforward  bit-map  processing  and  relatively 
small,  reliable  special  purpose  microprocessors  can  be  used,  to  support  these 
functions. 

IV.  CONCLUSIONS 

The  ASR/MTD  and  ARSR/MTD  radar  systems  can  be  used  for  weather  data 
extraction,  even  though  their  beam  patterns  and  scan  rates  are  tailored  for 
aircraft  surveillance  rather  than  radar  meteorology.  Based  on  our  studies 
using  clutter  from  BVA  (chapter  III)  and  the  MTD-II  experience  at  FAATC  and 
Burlington,  ,  we  have  concluded  that  ground  clutter  can  be  automat¬ 

ically  suppressed  without  "seriously"  compromising  the  weather  extraction 
function,  using  low-order  digital  filtering.  Spatial/temporal  smoothing  and 
contouring  algorithms  developed  for  the  MTD  two-level  weather  extraction  pro¬ 
cessor  will  support  the  six-level  processor  (five-thresholds).  We  see  no 
advantage  (and,  several  disadvantages)  in  combining  MTD  and  weather  extraction 
processing  in  a  single  unit,  and  recommend  that  they  be  implemented  separ¬ 
ately. 

The  key  concepts  involved  in  the  candidate  weather  processing  system  are: 

(1)  use  of  a  static  ground  clutter  map  to  select  the  form  of  clutter 
rejection  to  be  used  in  a  given  range-azimuth  cell 

(2)  iteration  of  the  ground  clutter  maps  in  a  variety  of  clear  weather 
situations  to  reduce  false  declarations 

(3)  assessment  of  clutter  rejection  techniques  based  on  actual  clutter 
data 

(4)  spatial/temporal  smoothing  and  contouring 

(5)  1/R2  STC  and  A/D  converter  of  I,  Q  sampling,  and 

(6)  an  appropriate  polarization  for  the  weather  channel. 

These  have  individually  been  demonstrated  over  the  past  few  years!  . 
However,  it  Dust  be  reemphasized  that  the  full  system  design  validation/ 
refinement  and  ATC  controller  evaluation  has  not  yet  been  carried  out. 
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It  should  also  be  noted  that  there  are  substantial  inaccuracies*  (e.g., 
due  to  weather  return  statistics,  breakdown  of  the  beam-filling  assumption 
and/or  bright  band  occurrence)  in  the  weather  return  estimates  which  one  can 
obtain  with  an  ASR  or  ARSR  even  when  clutter  is  not  present.  In  particular, 
the  operational  evaluation  should  consider: 

(1)  methods  of  identifying  and  flagging  situations  in  which  second  trip 
weather  echos  are  contained  in  the  estimated  first  trip  weather 
data,  and 

(2)  the  degree  to  which  the  use  of  a  larger  nmber  of  weather  threshold 
levels  for  the  ASR  or  ARSR  may  imply  a  greater  accuracy  in  the  esti¬ 
mates  than  is  in  fact  the  case. 


*Appendix  A  discusses  a  nunber  of  these  inaccuracies. 
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APPENDIX  A 


WEATHER  RADAR  PERFORMANCE 


1.  CANDIDATE  RADARS 

The  pertinent  performance  parameters  for  three  radars  of  Interest  are 
listed  In  Table  A-l.  NEXRAD  Is  a  hypothetical  next  generation  coherent  radar 
operating  at  S-Band.  It  may  be  thought  of  as  a  possible  contender  for  joint 
use  by  the  FAA,  the  Air  Force  Weather  Service  and  the  National  Weather  Service 
as  the  next  generation  en  route  weather  radar.  The  ARSRr A/MTD  (see  Fig.  A-l) 
should  be  viewed  as  the  present  ARSRr3  modified  by  the  type  of  front-end 
alterations  described  in  this  report  and  by  the  addition  of  MTD-type  process¬ 
ing.  Since  design  experience  and  field  data  has  been  obtained  on  a  radar 
similar  to  the  ASR-9/MTD,  performance  parasmters  have  been  included  for  it  as 
well. 

2.  PRECIPITATION  RETURN  COMPARISONS 

In  the  use  of  radar  for  precipitation  measurement  the  principal  proper¬ 
ties  have  been  the  total  returned  power  or  the  temporal  variation  of  that 
power  from  storm  regions.  It  is  generally  accepted  that  severe  storms  of  the 
type  that  present  significant  hasards  to  aviation,  contain  greater  quantities 
of  larger  water  droplets  and  thus  produce  larger  radar  returns.  The  received 
power  Pr  from  a  volume  of  precipitation  is: 

Pt  Gj  Gg  *2 

Pr  - -  ot 

(4*)3  R* 

ot  is  calculated  as  follows,  assuming  the  beam  to  be  filled  with  hydro¬ 
scatters  (“hydrometeors-) : 

(Reflective  scattering  volume)  x  (Scattering  coefficient,  n) 

Volume  ■  x  ( RB)  (R0)  .  x  .  Cx 

2  ~T  3Ln2  ~T 

for  Gaussian  Beams,  tw>-way,  and 

„  .  5-6  »  10-»  tl-6 

A‘ 

where: 

r  -  rain  rate  in  mm/hr. 
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COMPARISON  OF  Wx/No  FOR  CANDIDATE  RADARS 


+50  dB  +37dB  I  +34  dB 


BW  =  500  KHz 

LOSSES 

tR  2-  3  dB  RECEIVE  LINE  LOSSES 
Lt  »  2  dB  TRANSMIT  LINE  LOSSES 
Lp  *  1  dB  PROPAGATION  LOSS 

Fig.  A-l.  ARSR-4  front-end  performance  parameters  (assumed) . 


When  applying  this  equation  to  the  en  route  radar,  several  factors 
deserve  consideration.  The  first  is  the  statistical  nature  of  the  received 
power  from  precipitation.  A  single  measurement  is  likely  to  contain  large 
errors.  A  second  factor  is  the  beam-filling  assumption  used  in  the  derivation 
of  this  equation.  While  this  is  a  good  approximation  for  radars  with  a  narrow 
pencil  beam,  such  as  those  used  by  weather  radars,  for  a  fan  beam  radar  it 
may  underestimate  the  severity  of  the  storm  because  of  non-uniform  height 
profiles  or  low  altitude  storms  may  only  partially  fill  the  antenna  beam. 

This  problem  was  pointed  out  in  Kef.  1  by  Goonley  and  Hopson.  An  example 
of  storm  profile  data  is  shown  in  Fig.  A-2.  These  data  indicate  that  the 
average  power  from  integrating  over  the  coverage  height  will,  in  general, 
indicate  a  lower  z  value  than  the  maximum  actual  z  value  in  that  cell.  For 
example,  in  Figure  3  of  the  reference,  at  85  Km,  where  a  40  dBz  region  occurr¬ 
ed  between  4  to  6  Km  altitude,  the  volume-average  is  below  30  dBz. 

Another  factor  is  that  non-weather  related  targets  within  the  beam  cover¬ 
age  introduce  errors  in  the  measurement  accuracy.  For  the  en  route  radar  the 
need  to  provide  low  altitude  aircraft  coverage,  necessarily  introduces 
substantial  ground  clutter  return  energy  in  the  received  signal. 

The  processor  will  use  8  pulses  per  range  gate  to  estimate  Pr,  and  there¬ 
fore  z-level.  The  ARSR  will  operate  at  approximately  400  Hz,  thus  the  CPI 
interval  will  be  20  msec.  The  weather  signal  is  substantially  correlated 
pulse-to-pulse  (2.5  msec),  but  does  de-correlate  to  some  extent  over  the  CPI- 
interval.  The  extent  of  de-correlation  is  related  to  the  velocity  width  of 
the  return,  and  for  purposes  of  generating  a  typical  value  for  de-correlation, 
the  spectral  width  is  assumed  to  be  Gaussian  with  2  msec  variance,  yielding  an 
effective  doppler  width  of  43.5  Hz.  The  number  of  effective  independent 
samples  is  given  by 

Nef f  -  2WT 

where  W  •  43.5  Hz  and 
T  ■  20  msec 

Thus  there  are  approximately  two  Independent  samples  per  gate  per  CPI.  Eight 
range  gates  will  be  summed,  and  it  is  assumed  that  these  samples  are  essen¬ 
tially  uncorrelated,  yielding  the  equivalent  of  16  Independent  samples  for 
estimating  the  z-level.  Using  either  voltage  or  power  summation  (approxi¬ 
mately)  the  one-standard  deviation  of  the  estimate  will  be  about  1  dB .  Die 
estimate  variance  will  be  somewhat  greater  for  narrower  width  weather  returns, 
and  somewhat  smaller  for  wider  width  weather  returns.  Scair-to-scan  and 
spatial  smoothing  improve  the  estimator  to  some  extent,  and  our  experience 
with  the  ASR-MTD  two-level  processor  indicates  that  this  form  of  z-level 
intensity  estimation  will  be  adequate. 
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3.  CLUTTER/WEATHER 


Figures  A-3  and  A-4  provide  some  perspective  on  the  dynamic  range  of  the 
signals  Involved  and  the  relationship  between  ground  clutter  and  weather  re¬ 
turns  for  the  candidate  ARSR  and  NEXRAD  systems  likely  to  be  operated  at  L- 
band  and  S-band  by  the  FAA  in  the  future.  Figure  A-3  shows  parametric  curves 
of  S/N  for  weather  and  ground  clutter  returns  for  the  ARSR-4/MTD  operating  at 
L-Band  and  the  NEXRAD,  a  next  generation  coherent  weather  radar,  operating  at 
S-Band.  Fig.  A-4  shows  the  Clutter/Weather  ratios  for  these  systems  for  spe¬ 
cific  but  important  parametric  conditions.  The  weather  level  chosen  is  the 
upper  bound  of  the  NWS  weather  level  1  region,  and  would  be  the  lowest  level 
processed  for  extraction  by  the  ARSR  system.  Ihe  ground  clutter  back  scatter 
coefficient  of  -40  dB0  is  typical  of  much  terrain  that  would  be  visible  to 
en  route  weather  radars:  The  important  features  of  these  curves  are: 

a.  Very  wide  dynamic  range,  especially  when  -10  dB0  clutter  is 
considered,  and  weather  levels  greater  than  4-70  dBz. 

b.  With  the  beams  pointed  at  their  lowest  useful  elevation  angles 
(for  long  range  coverage)  Ci/Wx  ratios  are  high,  and  some 
clutter  mitigation  is  required,  in  order  to  extract  weather 
information  without  attendant  ground  clutter  false  alarms. 

c.  The  ARSR  problem  is  more  difficult  at  this  elevation  angle, 
and  it  stays  that  way,  whereas  the  NEXRAD  system  does  tilt  up 
and  gets  substantial  relief  in-close  from  typical  ground  clut¬ 
ter.  At  1.5°  elevation  (for  NEXRAD)  with  a  -25  dB  first  side- 
lobe  (one-way)  the  close-in  ground  clutter  problem  for  +30  dBz 
weather  is  non-existent.  The  ARSR  operates  at  this  elevation 
for  enhanced  long-range  A/C  coverage  and  must  cope  with  the 
ground  clutter  on  each  scan. 

Another  feature  shown  in  fig.  A-3  is  the  ARSR  Cl/No  after  typical 
1/R^-STC  front-end  attenuation  (typical  of  BVA-MTD  operation).  In  the 
20-100  Km  the  Ci/No  is  approximately  +30  dB.  Fig.  A-5  shows  a  +30  dB 
map  of  BVA,  obtained  using  a  dump  of  the  temporally  smoothed  ground  clutter 
map.  Ci/No  exceeds  +30  dB  within  the  closed  contours,  and  therefore  repre¬ 
sents  regions  exhibiting  higher  than  -40  dB0  ground  clutter  reflectivity. 
These  regions  therefore,  in  the  absence  of  weather,  exceed  the  30  dBz  weather 
level  threshold  by  at  least  18-25  dB  (it  would  by  10-17  dB  for  NEXRAD),  and 
zero-velocity  clutter  suppression  is  required  to  avoid  false  weather  display. 
For  this  particular  site  much  of  the  hilly  terrain  subtends  negative  elevation 
angles  to  the  radar,  and  NEXRAD  would  have  a  smaller  percentage  of  the  area 
exceeding  this  threshold  level.  Beyond  50  nmi  at  BVA,  very  few  range-azimuth 
cells  will  require  clutter  mitigation,  even  at  the  lowest  weather  threshold 
(30  dBz).  A  chart  of  the  NWS  weather  levels  is  shown  in  Fig.  A-6.  Fbr  the 
highest  weather  threshold  (57  dBz),  only  a  few  singularities  would  require 
clutter  filtering  or  possibly  fixed-map  censoring. 
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Fig.  A-3.  Signal-to-noise  vs  range  for  weather  and  ground  clutter. 
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A-6.  Categories  of  radar  Wx  echo  Intensity  (as  related  to  NWS  levels) 
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Another  feature  which  should  be  apparent  from  the  BVA  clutter  map  is  that 
most  of  the  features  are  small  in  area  compared  to  typical  weather  syteas,  and 
even  though  clutter  filtering  will  suppress  and  distort  some  weather  signals 
(near  zero-velocity) ,  spatial  continuity  considerations  should  allow  for  rela¬ 
tively  accurate  contouring  and  smoothing,  even  for  the  lowest  thresholded 
level. 


Present  NWS  weather  radar  systems  depend  on  intervention  by  trained 
weather  radar  meteorologists  to  identify  these  features  as  ground  clutter 
rather  than  rain,  but  the  ARSR-MTD  system  will  require  automated  suppression. 
The  NEKRAD  system  will  also  suppress  ground  clutter  automatically. 

4.  ARSR  OPERATION 


The  ARSR-()/MTD,  is  used  for  long  range  (200  nmi)  air  route  surveillance, 
and  uses  a  sector  beam  antenna  (1.5°  azimuth  by  5°  elevation  with  CSC* 
modification)  rotating  at  30°/sec.  (12  seconds  per  scan).  The  narrow  azimuth 
beam  and  high  scan  rate  broaden  the  ground  clutter  spectrum,  and  this  compli¬ 
cates  the  problem  of  clutter  mitigation,  since  weather  spectrum  which  overlaps 
must  be  suppressed  as  well.  The  large  beam  (compared  to  the  NEXRAD  1°  pencil 
beam)  will  generally  not  be  uniformly  filled  with  hydrometeors  (especially  at 
long  range)  and  the  estimate  of  rainfall  rate  will  be  low  as  a  consequence. 
However,  for  purposes  of  this  treatment,  beam-filling  is  assumed  at  all 
ranges . 


The  radar  will  be  operated  in  a  dual-staggered  PRF  mode  to  support  MTD 
processing,  with  each  coherent  processing  interval  (CPI)  consisting  of  at 
least  8  pulses.  Weather  will  be  processed  on  alternate  CPIs.  Considerations 
of  weather  level  estimate  variance  will  be  based  on  8  pulse  CPIs ,  as  presently 
tested  with  MTD  weather  extraction  processors.  The  CPIs  will  be  synchronized 
in  azimuth  to  support  MTD  clutter  map  maintenance,  and  will  permit  scan-to- 
scan  smoothing  in  range-azimuth  space.  The  weather  processor  will  be  synchro¬ 
nized  with  the  MTD  and  probably  use  the  same  range  gate  timing  of  1/8  nmi 
(1.54  u  sec)  sampling.  Thus,  although  the  separate  weather  extraction  pro¬ 
cessor  is  independent  of  the  MTD,  there  are  timing  constraints  imposed  by  MTD 
signalling/processing  strategy.  The  1.5°  elevation  angle  of  the  ARSR  antenna, 
used  in  this  discussion,  is  lower  than  present  systems  (normally  2.5°  or  high¬ 
er),  but  would  probably  be  used  with  MTD  processing,  since  clutter  mitigation 
can  be  managed,  and  this  gains  about  4  dB  in  sensitivity  for  long  range  and 
low  altitude  aircraft.  It  does,  of  course,  impose  an  additional  burden  on  the 
weather  extraction  processor. 
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MULTIPLE  TRIP  WEATHER  ECHOS 


One  addiCional  problem  common  Co  radar  measurement  methods  using  rela¬ 
tively  high  and  constant  PRF's  is  that  of  second-t ime-around  (STA)  echoes, 
arising  from  longer  range  targets  whose  round  trip  echo  time  exceeds  the  in¬ 
terpulse  period.  For  a  1200  PRF  characteristic  of  an  ASR,  this  unambiguous 
range  limit  is  approximately  68  nmi.  Thus,  possible  weather  returns  from 
distant  storms  in  the  90  to  110-nmi  range  interval,  for  example,  will  be 
superimposed  onto  desired  returns  between  20-  and  40-nmi  range.  The  resultant 
reflectivity  display  can  lead  to  highly  erroneous  conclusions  regarding  the 
actual  weather  situation.  Examples  of  this  on  an  ASR  in  Oklahoma  are  shown  in 
Zittel  [7]. 

Several  methods  have  been  suggested  for  mitigating  this  problem: 

(1)  comparison  of  the  weather  display  at  different  PRF's  [7], 

(2)  recognition  of  the  characteristic  elongated  shape  of  STA  echos  [7], 

(3)  use  of  a  low  PRF  (and  appropriate  clutter  rejection  strategy)  during 
periods  where  STA  echos  may  be  present  [8],  and 

(4)  decorrelation  of  second  trip  weather  returns  by  use  of  random  phase 
transmitter  changes  combined  with  Doppler  filtering  [4,7]. 

The  suggested  ASR/MTD  weather  channel  was  not  tested  in  the  experimental 
programs  at  FAATC,  BVD  and  BVT;  thus  STA  echos  were  not  a  significant  factor. 
However,  determining  an  appropriate  combination  of  hardware  features  and 
operational  use  methods  to  avoid  significant  weather  interpretation  errors 
should  be  a  key  element  of  the  ASR/MTD  "weather  channel”  validation  program. 
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APPENDIX  B 


CANDIDATE  CLUTTER  FILTERS 


The  class  of  filters  which  were  considered  for  the  weather  extraction 
processor  was  limited  to  those  which  could  be  realized  with  the  MTD  dual-PRF 
signaling  strategy.  In  all  cases,  eight  resulting  samples  were  to  be  avail¬ 
able  after  filtering,  to  support  the  z-level  estimators.  The  all-pass  and 
mean-leve 1-subtractor  (DC-removal  filter)  were  described  in  Section  III,  and 
are  essentially  trivial.  More  complex  FIR  filters  were  examined,  to  obtain 
more  or  less  clutter  mitigation,  without  "totally"  eliminating  the  weather 
signal.  Two  4-pulse  (with  feedback)  cancellors  are  shown  in  Figs.  B-l  and 
B-2,  and  standard  2-pulse  and  3-pulse  cancellors  are  shown  in  Figs.  B-3  and 
B-4.  The  ASR-MTD  high-pass  filter  response  is  shown  in  Fig.  B-5  for 
comparison,  but  it  is  not  a  candidate  for  this  processor. 

The  filter  set  used  for  the  example  described  in  Section  III  of  the  body 
of  this  report  was: 

Type  Cancellation 

0  -  0  dB 

1  —12  dB 

2  —20  dB 

3  —35  dB 

For  those  cells  requiring  greater  clutter  mitigation  than  is  possible 
with  these  candidate  filters,  it  is  recommended  that  fixed-map  censoring  be 
used,  as  was  necessary  for  some  extreme  cells  at  BTV,  in  the  MTD  processor. 


B-l 


Fig.  B-l.  4-pulse  equiripple  finite  impulse  response  filter. 
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Fig.  B-2.  4-pulse  equiripple  finite  impulse  response  filter. 
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Fig.  B-3.  Standard  2-pulse  canceller. 
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Fig.  B-4.  Standard  3-pulse  canceller. 
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APPENDIX  C 


LIST  OF  ACRONYMS  AND  ABBREVIATIONS 


A/C  Aircraft 

A/D  Analog-to-Digltal 

ARSK  Air  Route  Surveillance  Radar 

ARTCC  Air  Route  Traffic  Control  tenter 

ASR  Airport  Surveillance  Radar 

ATC  Air  Traffic  Control 

BPF  Ba.'d  Pass  Filter 

BVA  Bedford,  Virginia 

BW  Band  Width 

CFAR  Constant  False  Alarm  Rate 

Cl/No  Ratio  of  Clutter  return  in  the  i-th  cell  to  thermal  noise 

COHO  Coherent  Oscillator 

CP  Circular  Polarization 

CPI  Coherent  Processing  Interval 

dBsm  Decibels  with  respect  to  1.0  sq.  meter 

dBz  Decibels  with  respect  to  radar  reflectivity  factor,  z. 

DUP  Duplexer 

EL  Elevation  (angle) 

FAATC  Federal  Aviation  Administration  Technical  Center,  Atlantic 
City,  NJ 

FIR  Finite  Impulse  Response 

GMT  Greenwich  Mean  Time 

IF  Intermediate  Frequency 

I&Q  In-phase  and  Quadrature-phase 

MTD  Moving  Target  Detect(or) (ion) 

MTI  Moving  Target  Indicator  (ion) 

N EXRAD  Designates  a  hypothetical  "next  en-route  radar  design" 

NSSL  National  Severe  Storms  Laboratory 

NWS  National  Weather  Service 


APPENDIX  C  (CONT'D) 

LIST  OF  ACRONYMS  AND  ABBREVIATIONS  (CONTINUED) 


K/T  Receiver/ Transmit ter 

S/N  Signal-to-Noiae  Ratio 

STALO  Stabilized  Local  Oscillator 

STC  Sensitivity  Time  Control 

WFMU  Weather  and  Fixed  Map  Unit 

WG  Waveguide 

Mx  Weather 

Wx/Ci  Ratio  of  weather  return  to  clutter  return  in  the  i-th  cell 
ZVF  Zero  Velocity  Filter 


